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Nanoplasmonics has recently experienced explosive development with many novel ideas and dra- 
matic achievements in both fundamentals and applications. The spaser has been predicted and 
observed experimentally as an active element - generator of coherent local fields. Even greater 
progress will be achieved if the spaser could function as a ultrafast nanoamplifier - an optical coun- 
terpart of the MOSFET (metal-oxide-semiconductor field-effect transistor). A formidable problem 
with this is that the spaser has the inherent feedback causing quantum generation of nanolocalized 
surface plasmons and saturation and consequent elimination of the net gain, making it unsuitable 
for amplification. We have overcome this inherent problem and shown that the spaser can per- 
form functions of an ultrafast nanoamplifier in two modes: transient and bistable. On the basis 
of quantum density matrix (optical Bloch) equations we have shown that the spaser amplifies with 
gain > 50, the switching time < 100 fs (potentially, ~ 10 fs). This prospective spaser technology 
will further broaden both fundamental and applied horizons of nanoscience, in particular, enabling 
ultrafast microprocessors working at 10 — 100 THz clock speed. Other prospective applications are 
in ultrasensing, ultradense and ultrafast information storage, and biomedicine. The spasers are 
based on metals and, in contrast to semiconductors, are highly resistive to ionizing radiation, high 
temperatures, microwave radiation, and other adverse environments. 



I. INTRODUCTION 

Nanoplasmonic phenomena (see, e.g., Ref. Q]) unfold on 
the spatial scale between the skin depth l s w 25 nm (in 
noble metals) and the nonlocality radius £„; ~ vf ~ 1 
nm, where vp is the electron speed at the Fermi sur- 
face, and u> is optical frequency. Nanoplasmonics is ultra- 
fast: the temporal scale of the nanoplasmonic phenom- 
ena is between the coherent time of hundred attoseconds 
defined by the inverse bandwidth of the plasmonic fre- 
quency range (between uv and mid ir for plasmonic met- 
als and doped semiconductors )2£ and the surface plas- 
mon (SP) relaxation time 7" 1 ~ 10 — 100 fs (for no- 
ble metals in the visible to near-ir frequency range) 4 . 
Not just a promise anymore 5 , nanoplasmonics has deliv- 
ered a number of important applications: ultrasensing 6 , 
scanning near-field optical microscopy^, SP-enhanced 
photodetectors--, thermally assisted magnetic recording 9 , 
generation of extreme m*i°-, biomedical tests&ii, SP- 
assisted thermal cancer treatment 1 ^, and many others. 

To continue its vigorous development, nanoplasmonics 
needs an active device - near-field generator and ampli- 
fier of nanolocalized optical fields, which has until re- 
cently been absent. A nanoscale amplifier in microelec- 
tronics is the metal-oxide-semiconductor field effect tran- 
sistor (MOSFET^iii, which has enabled all contempo- 
rary digital electronics, including computers and com- 
munications and formed the present day technology as 
we know it. However, the MOSFET is limited by fre- 
quency and bandwidth to < 100 GHz, which is already a 
limiting factor in further technological development. An- 
other limitation of the MOSFET is its high sensitivity to 
temperature, electric fields, and ionizing radiation, which 
limits its use in extreme environmental conditions. 



An active element of nanoplasmonics is the SPASER 
(Surface Plasmon Amplification by Stimulated Emission 
of Radiation) that was proposed^ (see also Ref. [H) as a 
nanoscale quantum generator of nanolocalized coherent 
and intense optical fields. It has recently been observed 
experimentally 16 . However, there is a formidable prob- 
lem to set the spaser as a nanoscale quantum amplifier. 
A principal difference between a laser (quantum optical 
oscillator) and a quantum optical amplifier is in that the 
amplifier does not possess feedback: it does not have 
mirrors and any parasitic feedback due to light scatter- 
ing should be carefully minimized to prevent spontaneous 
generation. In a sharp contrast, in the spaser, as we dis- 
cuss below in this Introduction, the feedback is always 
inherently present because the metal plasmonic nanopar- 
ticle (spaser's core) supports SP modes whose fields ex- 
ert periodic perturbation on the gain medium causing 
the feedback, which principally cannot be removed. This 
feedback will cause the SP generation and the subsequent 
saturation of the gain. For any spaser (or any laser, for 
that matter) in a stationary (CW) regime the net ampli- 
fication must be zero (which is a condition of the stable 
CW operation). Therefore, one might assume that the 
spaser cannot serve as a nanoamplifier. 

In this article we have solved this problem suggesting 
two approaches to setting the spaser as a nanoamplifier: 
(i) Dynamic or transient approach, which based on the 
fact that during a femtosecond transient process after the 
population inversion is created but before the CW regime 
is established, the spaser possesses a net amplification, 
and (ii) Bistability approach that is based on the inclu- 
sion of a saturable absorber in the gain medium, which 
prevents spontaneous generation and sets the spaser as a 
bistable (logical) ultrafast nanoamplifier. This is done 
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employing a quantum-kinetic theory that is based on 
density-matrix equations, which are the adaptation of the 
optical Bloch equations of the laser theory for the case of 
the spaser. From these equations we have also obtained 
new and potentially useful results describing not only the 
ultrafast dynamics of the spaser but also its CW mode: 
"spasing curve" (that is a relation between the SP pop- 
ulation of the spasing mode and the pumping rate) and 
also the linewidth of the spaser. These results provide 
an excellent qualitative explanation of the recent experi- 
mental data^. 

On the basis of the present theory, one may envision 
femtosecond-cycle nanoplasmonic chips with a high de- 
gree of integration where spasers communicate and con- 
trol each other through their local optical fields or are 
connected with nanoplasmonic wires. These can perform 
ultrafast microprocessor functions. The spasers can also 
be integrated with nano-photodetectors and nanosensors 
to perform complex functions of intelligent ultrafast de- 
tection and sensing. In contrast to semiconductor tech- 
nology, the spasers are based on metals and, therefore, 
are highly resistive to ionizing radiation, high tempera- 
tures and other adverse environments, with possible ap- 
plications in space, nuclear industries and defense. These 
functions of spaser will further widen both the fundamen- 
tal and applied horizons of nanoplasmonics and, gener- 
ally, science and technology. 



A. Spaser Fundamentals 

We have recently discussed the physical principles of 
the spasingi 5 . but will briefly reiterate them here for the 
sake of completeness and self-containment. The spaser 
is a nanoplasmonic counterpart of the laser— The 
laser has two principal elements: resonator (or cavity) 
that supports photonic mode(s) and the gain (or active) 
medium that is population-inverted and supplies energy 
to the lasing mode(s). An inherent limitation of the laser 
is that the size of the laser cavity in the propagation di- 
rection is at least half wavelength and practically more 
than that even for the smallest lasers develope d 17 ' 18 ] 19 . 
In the spaser— this limitation is overcome. The spas- 
ing modes are surface plasmons (SPs) whose localiza- 
tion length is on the nanoscale 2 ^ and is only limited by 
the minimum inhomogeneity scale of the plasmonic metal 
and the nonlocality radius 21 l n \ ~ 1 nm. So, the spaser is 
truly nanoscopic, whose minimum total size can be just 
a few nm. 

The resonator of a spaser can be any plasmonic metal 
nanoparticle whose total size R is much less than the 
wavelength A and whose metal thickness is between l n [ 
and l s , which supports a SP mode with required fre- 
quency uj n . This metal nanoparticle should be sur- 
rounded by the gain medium that overlaps with the spas- 
ing SP eigenmode spatially and whose emission line over- 
laps with this eigenmode spectrally 4 . As an example, we 
consider a model of a nanoshell spase r 15 ' 22 , which is illus- 
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FIG. 1: Schematic of spaser geometry, local fields, and fun- 
damental processes leading to spasing. (a) Nanoshell geom- 
etry and the local optical field distribution for one SP in an 
axially-symmetric dipole mode. The nanoshell has aspect ra- 
tio r\ — 0.95. The local field magnitude is color-coded by the 
scale bar in the right-hand side of the panel (b) The same as 
(a) but for a quadrupole mode, (c) Schematic of a nanoshell 
spaser where the gain medium is outside of the shell, on the 
background of the dipole-mode field, (d) The same as (c) but 
for the gain medium inside the shell, (e) Schematic of the 
spasing process. The gain medium is excited and population- 
inverted by an external source, as depicted by the black arrow, 
which produces electron-hole pairs in it. These pairs relax, as 
shown by the green arrow, to form the excitons. The exci- 
tons undergo decay in the ground state emitting SPs in the 
nanoshell. The plasmonic oscillations of the nanoshell stim- 
ulates this emission, supplying the feedback for the spaser 
action. 



trated in Fig. [T] Panel (a) shows a silver nanoshell car- 
rying a single SP (plasmon population number N n = 1) 
in the dipole eigenmode. It is characterized by a uniform 
field inside the core and hot spots at the poles outside 
the shell with the maximum field reaching ~ 10 6 V/cm. 
Similarly, Fig. [1] (b) shows the quadrupole mode in the 
same nanoshell. In this case, the mode electric field is 
non-uniform, exhibiting hot spots of ~ 1.5 x 10 6 V/cm of 
the modal electric field at the poles. These high values of 
the modal fields is the underlying physical reason for very 
strong feedback in the spaser. Under our conditions, the 
electromagnetic retardation within the spaser volume can 
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be safely neglected. Also, the radiation of such a spaser is 
a weak effect: the decay rate of plasmonic eigenmodes is 
dominated by the internal loss in the metal. Therefore, it 
is sufficient to consider only quasistatic eigenmode o 20 ' 23 
and not their full electrodynamic counterparts 24 , which 
also is necessary for the ultrafast operation of the spaser 
as a nanoamplifier. 

Note that for high aspect ratios (say, 77 = 0.95 consid- 
ered in some of the examples) and relatively small radii 
(R = 12 nm) the thickness of the metal shell may be 
rather small (less than 1 nm), which may lead to nonlo- 
cal effects including Landau damping. In this case, the 
radius of the spaser can be safely increased by a factor 
of 2 to 3 (this will not adversely effect the spasing as the 
spaser is fully scalable, as we show below in Sec. Ill Bl) . 
Also one should keep in mind that the nonlocal effects 
in thin nanoshclls arc small due to the special symme- 
try reasons^ 5 -. We also consider a number of examples 
with significantly higher aspect ratios where the metal 
thickness is a few nanometers or larger. 

For the sake of numerical illustrations of our theory, 
we will use the dipole eigenmode [Fig. Q] (a)]. There are 
two basic ways to place the gain medium: (i) outside the 
nanoshell, as shown in panel (c), and (ii) in the core, as in 
panel (d) , which was originally proposed in Ref . |22| . As 
we have verified, these two designs lead to comparable 
characteristics of the spaser. However, the placement of 
the gain medium inside the core illustrated in Fig. [1] (d) 
has a significant advantage because the hot spots of the 
local field are not covered by the gain medium and are 
sterically available for applications. 

Note that any Z-multipole mode of a spherical particle 
is, indeed, 21 + 1-times degenerate. This may make the 
spasing mode to be polarization unstable, like in lasers 
without polarizing elements. In reality, the polarization 
may be pinned and become stable due to deviations from 
the perfect spherical symmetry, which exist naturally or 
can be introduced deliberately. More practical shape for 
a spaser may be a nanorod, which has a mode with the 
stable polarization along the major axis. However, a 
nanorod is a more complicated geometry for theoretical 
treatment, and we will consider it elsewhere. 

The level diagram of the spaser gain medium and 
the plasmonic metal nanoparticle is displayed in Fig. 
Q] (e) along with a schematic of the relevant energy 
transitions in the system. The gain medium chro- 
mophores may be semiconductor nanocrystal a 4 ' 26 , dye 
molecule s 27 ' 28 , rare-earth ions^ 2 -, or electron-hole excita- 
tions of a bulk semiconducto r 17 ! 18 . For certainty, we will 
use a semiconductor-science language of electrons and 
holes. 

The pump excites electron-hole pairs in the chro- 
mophores [Fig. [1] (e)], as indicated by the vertical black 
arrow, which relax to form excitons. The excitons con- 
stitute the two-level systems that are the donors of en- 
ergy for the SP emission into the spasing mode. In vac- 
uum, the excitons would recombine emitting photons. 
However, in the spaser geometry, the photoemission is 



strongly quenched due to the resonance energy trans- 
fer to the SP modes, as indicated by the red arrows in 
the panel. The plasmons in the spaser mode create the 
high local fields that excite the gain medium and stim- 
ulate more emission to this mode, which is the feedback 
mechanism. If this feedback is strong enough and the 
life time of the spaser SP mode is long enough, then an 
instability develops leading to the avalanche of the SP 
emission in the spasing mode and spontaneous symmetry 
breaking, establishing the phase coherence of the spasing 
state. Thus the establishment of spasing can be called 
a nonequilibrium phase transition, as in the physics of 
lasing. 



B. Brief Overview of Latest Progress in Spasers 
and Nanolasers 

After the original theoretical proposal and predic- 
tion of the spaser—, there has been an active develop- 
ment in this field, both theoretical and experimental. 
We comment below only on some representative pub- 
lications. Among theoretical developments, a nanolcns 
spaser has been proposed 2 ^,, which possesses a nanofocus 
("the hottest spot") of the local fields. In Refs.l4ll29[ only 
the necessary condition of spasing has been established 
on the basis of the perturbation theory. 

There have been theories published describing the 
spaser (or, "nanolaser" as sometimes it is called) phe- 
nomenologically, on the basis of classic linear elec- 
trodynamics by considering the gain medium as a 
dielectric with a negative imaginary part of the 
permittivit y 22 ' 30 ! 31 ' 32 . Such electrodynamic approaches 
do not take into account the nature of the spasing as a 
spontaneous symmetry breaking. This leads to principal 
differences of their results from the present microscopic 
quantum-mechanical theory in the region of spasing, as 
we discuss below in Sec. IIII Al in conjunction with Fig. 

Note that before the spasing threshold, the above- 
mentioned phenomenological theories are, in principle, 
applicable. There has also been a theoretical publication 
on a bowtie spaser (nanolaser) with electrical pumping 33 . 
It is based on balance equations and only the CW spasing 
generation intensity is described. Yet another theoretical 
development has been a proposal of the lasing spaser—, 
which is made of a plane array of spasers. The theoreti- 
cal publications mentioned above in this paragraph deal 
with the CW regime. In contrast, in this article we are 
most interested in the ultrafast behavior of the spaser. 

There have also been a theoretical proposal of a spaser 
( "nanolaser" ) consisting of a metal nanoparticle coupled 
to a single chromophoro 3 ^. In this paper, a dipole-dipole 
interaction is illegitimately used at very small distances 
r where it has a singularity (diverging for r — > 0), lead- 
ing to a dramatically overestimated coupling with the 
SP mode. As a result, a completely unphysical predic- 
tion of CW spasing due to single chromophore has been 
obtained 35 . In contrast, our theory is based on the full 
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(exact) field of the spasing SP mode. As our results of 
Sec. lIIII below show, hundreds of chromophores per metal 
nanoparticle are realistically requited for the spasing even 
under the most favorable conditions. 

There has been a vigorous experimental quest toward 
the spaser. Stimulated emission of SPPs has been ob- 
served in a proof-of-principle experiment using pumped 
dye molecules as an active (gain) medium^!. There have 
also been later experiments that demonstrated strong 
stimulated emission compensating a significant part of 
the SPP los o 28 ! 36 ^^. As a step toward the lasing 
spaser, the first experimental demonstration has been 
reported of compensating Joule losses in a metallic pho- 
tonic metamaterial using optically pumped PbS semicon- 
ductor quantum dots2£. In another development, a tun- 
able free-electron light nanosource ("light well") has been 
demonstrated^, which is based on spontaneous emis- 
sion in a periodically layered metal-dielectric structure. 
There have also been experimental investigations report- 
ing the stimulated emission effects of SPs in plasmonic 
metal nanoparticles surrounded by gain media with dye 
molecule a 40 i 41 . An electrically-pumped nanolaser with 
semiconductor gain medium has been demonstrated 18 
where the lasing modes are SPPs. A nanolaser with 
optically-pumped semiconductor gain medium and a 
hybrid semiconductor-silver SPP waveguide has been 
demonstrated with an extremely tight transverse mode 
confinement 1 ^. 

Finally, an observation has been reported of a true 
spaser by M. Noginov et al^. This spaser is a chemically 
synthesized gold nanosphere of radius 7 nm surrounded 
by a dielectric shell of a 21 nm outer radius containing im- 
mobilized dye molecules, which is sometimes referred to 
as "Cornell nanodot" . Under nanosecond optical pump- 
ing in the absorption band of the dye, this spaser de- 
velops a relatively narrow-spectrum and intense visible 
emission that exhibits a pronounced threshold in pump- 
ing intensity. The observed characteristics of this spaser 
are in an excellent qualitative agreement and can be fully 
understood on the basis of the corresponding theoretical 
results obtained below in Sec. IIII Al 

This article is organized as the following. In Sec. Ill Al 
for the case of the spaser, we formulate quantum den- 
sity matrix equations (also called in the literature the 
optical Bloch equations). In Sec. Ill Bl we consider these 
equations for a stationary (CW) case where we find the 
conditions of spasing, the population of the SPs in the 
spasing mode, and the spasing line shape. In Sec. IIII Al we 
illustrate the spasing curve and spectral composition of 
the spaser generation. In Sec. IIII Al we describe the CW 
spaser as a bistable nanodevice. The major results of 
this articles regarding the ultrafast kinetics of the spaser 
as a generator and quantum nanoamplifier are presented 
in Sec. [IVJ We briefly discuss the obtained results and 
conclude in Sec. [V] 



II. EQUATIONS OF SPASER 

A. Quantum Density Matrix (Optical Bloch) 
Equations for Spaser 

The SP eigenmodes ^ n (r) are described by a wave 
equation (with homogeneous boundary conditions) 4 -^ 

V6(r)V^„(r) = * n VVn(r) , (1) 

where n is the mode number, s n is corresponding eigen- 
value, and O(r) is the characteristic function equal to 
1 for r in the metal component and in the dielectric. 
Note that the eigenvalues s n are all real and contained 
in the range 1 > s n > 0. The eigenmodes are normal- 
ized by an integral over the volume V of the system, 
J v | Vy„(r)| 2 d 3 r = 1. The physical frequency uj n of the 
SPs is defined by an equation Re [s(w n )] = s n , where 
s(lj) — Sd/ [s<i — £m{u)] is Bergman's spectral parame- 
ter, Ed is the permittivity of the ambient dielectric, and 
e m (oj) is the metal permittivity. 

The electric field operator— of the quantized SPs is 4 - 

n \ £ dS n J 

(2) 

where and a n are the SP creation and annihilation 
operators, and s' n = Re [ds(ui n )/dui n ]- 
The spaser Hamiltonian has the form 

H = H g + hJ2^naia n - Yl E ( r p) d(p) . ( 3 ) 

n p 

where H g is the Hamiltonian of the gain medium, p is 
an index (label) of a gain medium chromophorc, r p is 
its coordinate vector, and d' p ' is its dipole moment op- 
erator. In this paper, we will treat the active medium 
quantum mechanically but the SPs quasiclassically, con- 
sidering classical quantity (c-number) a n with 
time dependence as a n — ao n exp(— iu)t), where ao n is a 
slowly-varying amplitude. The number of coherent SPs 
per spasing mode is then given by N p — |ao n | 2 . This 
approximation neglects the quantum fluctuations of the 
SP amplitudes. However, when necessary, we will take 
into account these quantum fluctuations, in particular, 
to describe the spectrum of the spaser. 

Introducing p^ as the density matrix of a pth chro- 
mophore, we can find its equation of motion in a conven- 
tional way by commutating it with the Hamiltonian §3§ as 
ihp^ = [p^ p \ H], where the dot denotes temporal deriva- 
tive. We will use the standard rotating wave approxima- 
tion (RWA), which only takes into account the resonant 
interaction between the optical field and chromophores. 
We denote |1) and |2) as the ground and excited states 
of a chromophore, with the transition |2) ^ |1) reso- 
nant to the spasing plasmon mode n. In this approx- 
imation, the time dependence of the nondiagonal ele- 
ments of the density matrix is (p^) 12 = P%2 exp(iu>t), 
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and (p^) 21 — P12 exp(— iwt), where pf 2 is a time- 
independent amplitude denning the coherence (polariza- 
tion) for the 1 2) ^ |1) spasing transition in a pth chro- 
mophore of the gain medium. 

Introducing a rate constant Ti2 to describe the polar- 



ization relaxation and a difference 



» 



population inversion on this spasing transition, we derive 
an equation of motion for the non-diagonal element of the 
density matrix as 



-<j>) 

Pl2 



[i (uj - LJ12) + 



ri 2 ]M ; 



(p) 
2 



*21 i2 12 



(4) 



where = —A n df 2 Vip n (r p )aQ n /h is the Rabi fre- 

quency for the spasing transition in a pth chromophore, 
and dj^ is the corresponding transitional dipole element. 

Note that always d^ is either real or can be made real 
by a proper choice of the quantum state phases, making 
the Rabi frequency Oj^ a iso a real quantity. 

An equation of motion for n 21 can be found in a stan- 
dard way by commutating it with H. To provide con- 
ditions for the population inversion (n 21 > 0), we imply 
existence of a third level. For simplicity, we assume that 
it very rapidly decays into the excited state |2) of the 
chromophore, so its own populations is negligible. It is 
pumped by an external source from the ground state (op- 
tically or electrically) with some rate that we will denote 
g. In this way, we obtain the following equation of mo- 
tion: 



i(p)i 



'21 



-41m 



P12 "21 



72 



1 



'21 



9 



1 



>) 
'21 



( 5 ) 

where 72 is the decay rate |2) — * |1). 

The stimulated emission of the SPs is described as 
their excitation by the coherent polarization of the gain 
medium. The corresponding equation of motion can be 
obtained using Hamiltonian ([3]) and adding the SP relax- 
ation with a rate of 7„ as 



d „ = [i(u- oj„ ) - In] a 0n + i V" p$*Sl^ 



p 



(6) 



Another relevant process is spontaneous emission of 
SPs by a chromophore into a spasing SP mode. The 
corresponding rate 7^ f° r a chromophore at a point r p 
can be found in a standard way using the quantized field 
© as 

(p) n A\ 2 (ri2+7n) 2 

ft 7n (Wl2 -Wn) + ^12 + In) 

(7) 

As in Schawlow- Towns theory of laser-line width, this 
spontaneous emission of SPs leads to the diffusion of the 
phase of the spasing state. This defines width 7 S of the 
spasing line as 



2(2N p + 1) 



(p) 

72 



(8) 



This width is small for a case of developed spasing when 
N p ^> 1. However, for N p ~ 1, the predicted width 
may be too high because the spectral diffusion the- 



ory assumes that 7 S 



< 



In- 



To take into account this 



limitation in a simple way, we will interpolate to find 
the resulting spectral width T s of the spasing line as 

r, = (7 I 7 2 + 7 s - 2 )" 1/2 . 

We will also examine the spaser as a bistable (logical) 
amplifier. One of the ways to set the spaser in such a 
mode is to add a saturable absorber. This is described by 
the same Eqs. (|4|)-([6]) where the chromophores belonging 
to the absorber are not pumped by the external source 
directly, i.e., for them in Eq. © one has to set 5 = 0. 

Numerical examples are given for a silver nanoshcll 
where the core and the external dielectric have the same 
permittivity of = 2; the permittivity of silver is 
adopted from Ref. |42|. The following realistic parameters 
of the gain medium are used (unless indicated otherwise) : 
d 12 = 1.5 x 1CT 17 esu, hT 12 = 10 meV, 72 = 4 x 10 12 s" 1 
(this value takes into account the spontaneous decay into 
SPs), and density of the gain medium chromophores is 
p = 2.4 x 10 20 cm -3 , which is realistic for dye molecules 
but may be somewhat high for semiconductor quantum 
dots that were proposed as the chromophores^ and used 
in experiments^. We will assume a dipole SP mode and 
chromophores situated in the core of the nanoshcll as 
shown in Fig. Q] (d). This configuration are of advantage 
both functionally (because the region of the high local 
fields outside the shell is accessible for various applica- 
tions) and computationally (the uniformity of the modal 
fields makes the summation of the chromophores trivial, 
thus greatly facilitating numerical procedures). 



B. Equations for Stationary (CW) Regime 

Physically, the spaser action is a result of spontaneous 
symmetry breaking when the phase of the coherent SP 
field is established from the spontaneous noise. Math- 
ematically, the spaser is described by homogeneous dif- 
ferential Eqs. (0J- ([6]). These equations become homo- 
geneous algebraic equations for the stationary [or con- 
tinuous wave (CW)] case. They always have a trivial, 
zero solution. However, when their determinant vanishes, 
they also possess a nontrivial solution describing spasing, 
whose condition is 



(0J S - U>21 + i^l2Y 



(9) 



p 



Ap) 
'-21 



(p) 
12 



A n d^Vip n (r p )/h is the single-plasmon 



where 

Rabi frequency. The population inversion of a pth chro- 
mophore 72,21 is explicitly expressed as 



» 



(.9 - 72) x 



(10) 



5 + 7 2 + 4iV n 



nip) 



(lu s -lu 2 i) +T 2 2 
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From the imaginary part of Eq. (|10p we immediately 
find the spasing frequency 



uj s = (7„w2i + ri 2 w„) / (7„ + r i2 ) 



(ii) 



which generally does not coincide with either the gain 
transition frequency W21 or the SP frequency u> n , but is 
between them (this is a frequency walk-off phenomenon 
similar to that of laser physics). Substituting Eq. (TTTj) 
back to Eqs. (p^|) - (|TT]) . we obtain a system of equations 



(in 



T 



12) 



7n r i2 [(^21 - LUnf + (^12 + InY 

n< 2i = i • 

n 21 = (.9~72) X 



(12) 



5 + 72 + 



f2 



(p) 



(r i2 + 7n) 



(wi2-o; n ) +(r i2 +7„) 



(13) 



This system defines the stationary (CW) number of SPs 
per spasing mode N n . 

Since nf^ < 1, from Eqs. (TT2")) . (|13p we immediately 
obtain a necessary condition of the existence of spasing, 



(7« + r 12 r 



7«ri 



(w 2 i 



2 + (r 12 + 7 „) 1 



El" 



> 1 . 



(14) 

This expression is fully consistent with Ref. |4|. The fol- 
lowing order of magnitude estimate of this spasing condi- 
tion has a transparent physical meaning and is of heuris- 
tic value, 



d\ 2 QN c 

nr 12 v n 



> i 



(15) 



where Q = u)/~f n is the quality factor of SPs, V n is the 
volume of the spasing SP mode, and N c is the of number 
of the gain medium chromophorcs within this volume. 
Deriving this estimate, we have neglected the detuning, 
i.e., set u)2i — w n = 0. We also used the definitions of A n 
of Eq. ^) and f2j^ given after Eq. (jTOJ) , and the estimate 
|V<yj„(r)| 2 ~ 1/V following from the normalization of the 
SP eigenmodes / |V(^„(r)| 2 d 3 r = 1 of Ref.[H The result 
of Eq. (fTS"|) is, indeed, in agreement with Ref. where it 
was obtained in slightly different notations. 

It follows from Eq. (fT5")) that for the existence of spas- 
ing it is beneficial to have a high quality factor Q, a 
high density of the chromophores, and a large transition 
dipole (oscillator strength) of the chromophore transi- 
tion. The small modal volume V n (at a given number 
of the chromophores iV c ) is beneficial for this spasing 
condition: physically, it implies strong feedback in the 
spaser. Note that for the given density of the chro- 
mophores p c — N c /V n , this spasing condition does not 



explicitly depend on the spaser size, which opens up a 
possibility of spasers of a very small size limited from the 
bottom by only the nonlocality radius l n i ~ 1 nm. An- 
other important property of Eq. (|15p is that it implies the 
quantum-mechanical nature of spasing and spaser ampli- 
fication: this condition essentially contains the Planck 
constant h and, thus, does not have a classical counter- 
part. Note that in contrast to lasers, the spaser theory 
and Eq. (| 1 5|) in particular do not contain speed of light, 
i.e., they are quasistatic. 



III. SPASER IN MODE OF CONTINUOUS 
WAVE NANOSCALE QUANTUM GENERATOR 

A. Kinetics of CW Spaser 

The "spasing curve", i.e., the dependence of the coher- 
ent SP population N n on the excitation rate g, obtained 
by solving Eqs. (fl2)) . (|13p . is shown in Fig. [2] (a) for four 
types of the silver nanoshells with the frequencies of the 
spasing dipole modes as indicated, which are in the range 
from near ir (tku s = 1.2 eV) to mid visible (hjj s = 2.2 eV). 
In all cases, there is a pronounced threshold of the spas- 
ing at an excitation rate gth ~ 10 12 s . Soon after the 
threshold, the dependence N n (g) becomes linear, which 
means that every quantum of excitation added to the ac- 
tive medium with a high probability is stimulated to be 
emitted as a SP, adding to the coherent SP population. 

While this is similar to conventional lasers, there is a 
dramatic difference for the spaser. In the lasers a similar 
relative rate of the stimulated emission is achieved at a 
photon population of ~ 10 18 — 10 20 , while in the spaser 
the SP population is N n < 100. This is due to the much 
stronger feedback in spasers because of the much smaller 
modal volume V n - see discussion of Eq. (fTS")) . The shape 
of the spasing curves of Fig. [5] (a) (the well-pronounced 
threshold with the linear dependence almost immediately 
above the threshold) is in an excellent qualitative agree- 
ment with the experiment^. Note that the recent SPP 
nanolaser with a tight transverse confinement of the las- 
ing modes does not exhibit a pronounced threshold^, in 
contrast to the spaser. 

The population inversion number n 2 i as a function of 
the excitation rate g is displayed in Fig. [2] (b) for the 
same set of frequencies (and with the same color cod- 
ing) as in panel (a). Before the spasing threshold, n 2 i 
increases with g to become positive with the onset of 
the population inversion just before the spasing thresh- 
old. For higher g, after the spasing threshold is exceeded, 
the inversion n 2 i becomes constant (the inversion pin- 
ning). The pinned levels of the inversion are very low, 
n 2 i ~ 0.01, which again is due to the very strong feed- 
back in the spaser. 

The spectral width T s of the spaser generation is due 
to the phase diffusion of the quantum SP state caused 
by the noise of the spontaneous emission of the SPs into 
the spasing mode, as described by Eq. ([5]). This width 
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FIG. 2: Spaser SP population and spectral characteristics 
in the stationary state. The computations are done for a 
silver nanoshell with the external radius R2 = 12 nm; the 
detuning of the gain medium from the spasing SP mode is 
H (u>2i — Un) = —0.02 eV. The other parameters are indicated 
in Sec. [TTJ (a) Number N n of plasmons per spasing mode 
as a function of the excitation rate g (per one chromophore 
of the gain medium). Computations are done for the dipole 
eigenmode with the spasing frequencies u s as indicated, which 
were chosen by the corresponding adjustment of the nanoshell 
aspect ratio, (b) population inversion ri\i as a function of the 
pumping rate g. The color coding of the lines is the same 
as in panel (a), (c) The spectral width T s of the spasing line 
(expressed as fiP s in meV) as a function of the pumping rate g. 
The color coding of the lines is the same as in panel (a), (d)- 
(f ) Spectra of the spaser for the pumping rates g expressed in 
the units of the threshold rate g t h, as indicated in the panels. 
The curves are color coded and scaled as indicated. 



is displayed in Fig. [2] (c) as a function of the pumping 
rate g. At the threshold, T s is that of the SP line j n 
but for stronger pumping, as the SPs accumulate in the 
spasing mode, it decreases oc N~ , as given by Eq. ([8]). 
This decrease of T s reflects the higher coherence of the 
spasing state with the increased number of SP quanta 
and, correspondingly, lower quantum fluctuations. This 
is similar to the lasers as described by the Schawlow- 
Townes theory^. 

The developed spasing in a dipole SP mode will show 
itself in the far field as an anomalously narrow and in- 
tense radiation line. The shape and intensity of this line 
in relation to the lines of the spontaneous fluorescence of 
the isolated gain medium and its SP-enhanced fluores- 
cence line in the spaser is illustrated in Figs. (d)-(f). 
Note that for the system under consideration, there is a 
20 meV red shift of the gain medium fluorescence with 



respect to the SP line center. It is chosen so to illustrate 
the spectral walk-off of the spaser line. For one percent 
in the excitation rate above the threshold of the spas- 
ing [panel (d)], a broad spasing line (red color) appears 
comparable in intensity to the SP-enhanced spontaneous 
fluorescence line (blue color). The width of this spasing 
line is approximately the same as of the fluorescence, but 
its position is shifted appreciably (spectral walk-off) to- 
ward the isolated gain medium line (green color). For 
the pumping twice more intense [panel (e)], the spaser- 
line radiation dominates, but its width is still close to 
that of the SP line due to significant quantum fluctua- 
tions of the spasing state phase. Only when the pumping 
rate is an order of magnitude above the threshold, the 
spaser line strongly narrows [panel (f)], and it also com- 
pletely dominates the spectrum of the radiation. This is 
a regime of small quantum fluctuations, which is desired 
in applications. 

These results in the spasing region are different in 
the most dramatic way from previous phenomcnolog- 
ical models, which are based on the consideration of 
a gain medium that has negative imaginary part of 
its permittivity plus lossy metal nanosystem, described 
purely electrodynamically 22 ^. For instance, in a "toy 
model"—, the width of the resonance line tends to zero 
at the threshold of spasing and then broadens up again. 
This distinction of the present theory is due the na- 
ture of the spasing as a spontaneous symmetry breaking 
(noncquilibrium phase transition) leading to the estab- 
lishment of the coherent SP state. This state is influ- 
enced by the phase relaxation due to the spontaneous 
SP emission into the spasing mode leading to the finite 
spectral width of the spasing line T s cx N~ x . There- 
fore the present theory requires a quantum mechanical 
consideration of the gain medium. Note that below the 
spasing threshold, the phenomenological theorie a 22 ' 32 are 
applicable, describing such important effects as loss com- 
pensation in metamaterials 44 . 



B. Bistable CW Spaser 

Bistables in microelectronics are MOSFET-based de- 
vices that have two stable states describing logical 1 and 
0. Their output "logical" level changes when the input 
exceeds a certain threshold value. Below we show that 
the spaser can operate as a bistable based on the quan- 
tum amplification. 

Generally, bistability is a result of nonlinearity in the 
system. We examine below the bistability in the spaser 
where the nonlinearity is due to the presence of a sat- 
urable absorber. Such an absorber is a chromophore 
whose absorption overlaps with the spasing line, but 
which does not directly absorb the radiation that pumps 
the spaser. Formally, it is described by Eqs. (fT2|) . (fT3]) 
where the pumping rate 5 = for the index p correspond- 
ing to the saturable absorber. We will assume that this 
absorber is distributed in the space the same way as the 
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FIG. 3: Bistability in spaser with saturable absorber in a sta- 
tionary spasing state, (a) Dependence of the SP population 
number N n in the spasing mode on the pumping rate g for dif- 
ferent concentrations p a of the saturable absorber. The curves 
are color coded corresponding to p a shown in the units of the 
concentration p of the active medium chromophores. The 
black line shows the threshold curve (separatrix) between the 
bistable and uniquely stable solutions, (b) The dependence 
of the population inversion ri2i on the pumping rate g for the 
spasing states. The color coding is the same as in panel (a). 
Note that the lower branches of the curves in this panel corre- 
spond to the upper ones in panel (a) and vice versa, (c) Same 
as in panel (a) but only the stable branches of the curves 
are displayed, illustrating the physical behavior of the spaser 
population N n . Note that one of the stable branches coin- 
cides with the horizontal axis. The vertical arrows illustrate 
the hysteretic transitions between the two stable branches for 
p a = 3p. (d) Same as in panel (c) but for the population in- 
version 7i2i. The black curve starting with ri2i = — 1 shows a 
trivial stable solution corresponding to the absence of spasing 
(N n = 0). 



gain chromophores but with a different density p a . 

Results of a numerical solution of Eqs. (fT2| . (fT3j) for 
different values of p a relatively to the concentration p 
of the gain medium chromophores are shown in Fig. [3J 
Panel (a) displays dependence of the SP population num- 
ber N n in the spasing mode on the pumping rate g, and 
the panel (b) shows the corresponding dependence of 
the population inversion 77,21 of the gain medium. We 
note first that there always exists also a trivial, non- 
spasing solution N n = 0. For values of g above critical 
(depending on p a ), there are also nontrivial solutions. 
For p a > 10 _3 p, these nontrivial solutions consist of two 
branches: high-iV n and low-AT n . Note that the high-iV n 
branch in panel (a) corresponds to the low-77,21 branch in 
panel (b), and vice versa. 

Thus it appears that there is a tri-stability. However, 
this impression is incorrect. The lower branches in Fig. 
131 (a) [and, correspondingly, the upper branches in panel 
(b)] describe unstable solutions that are not realizable 
physically. This can be understood already from the fact 
that along these branches the SP population decreases 
with increasing pumping, which is completely unphysical. 

It is a bistability, not the tri-stability, which takes place 



in actuality. This is illustrated in Figs. [3J (c) and (d) ob- 
tained by isolating the stable branches of the solutions. 
As one can see, with an increase of the pumping rate g, 
the solutions appear at critical pumping rates that in- 
crease with the saturable absorber concentration p a . As 
the critical g for a given p a is reached, the nonzero-iV n 
branch appears with a discontinuity. Both the N n = 
and N n > branches are stable and can keep their states 
indefinitely. The transition between these two stable 
states can be induced by either adding SP quanta to or 
removing them from the spasing mode, as illustrated by 
arrows for p a = 'ip. This shows that a spaser with a sat- 
urable absorber can serve functions of both a nanoscopic 
memory cell and a bistable based on quantum amplifica- 
tion. The dynamics of such a device is femtosecond as 
we show below in Sec. IIVI 



In the bistable case there is no perfect pinning of the 
population inversion 77,21, as Fig. [3J(d) shows. When the 
pumping rate g is increased, the system moves along the 
stable no-spasing (N n = 0) solution denoted by the black 
line, where the inversion 77,21 significantly overshoots its 
values for the spasing branches (colored lines). When the 
transition to spasing occurs, induced, e.g., by an injection 
of SPs quanta, its is always discontinuous, as the vertical 
arrows indicate. Along the spasing (N n > 0) branches 
(colored lines), the inversion counterintuitevely decreases 
with increased pumping due to the stimulated emission. 



It is instructive to compare the behavior of the spaser 
shown above in Fig. [3J (see also its discussion) with the 
corresponding behavior of lasers. As we have already 
mentioned in Sec. U there are similarities and also sub- 
stantial differences. A principal source of these differ- 
ences is due to the fact the spaser SP modes are much 
stronger localized than the photonic modes in conven- 
tional lasers and even the SPP modes in the nanolasers 
of Refs. Il8lll9l . The much smaller modal volume V n of the 
spasing SP modes causes much higher local fields of these 
SPs [cf. Fig. m and, consequently, much stronger feed- 
back in the spaser - cf. Eq. (fT5|) . Therefore the spaser 
operates, including the bistable operation, at many or- 
ders of magnitude smaller population of SPs (typically, 
N n ~ 100) per spasing mode as compared to the num- 
ber of quanta in the resonators of the conventional lasers 
(~ 10 18 — 10 20 ) and even the polaritonic nanolasers. Re- 
lated to this fact are the very high values of the local 
fields in the spaser ~ 10 6 ^/N n ~ 10 7 V/cm, as is pointed 
out above in this paragraph. These strong local fields and 
the strong feedback caused by them are also responsible 
for the pronounced linear dependence of the SP popu- 
lation N n on the pumping rate [see Fig. [3J (c)] and the 
pronounced pinning of the population inversion 77,21 of 
the gain medium [see Fig. [2] (b) and Fig. [3J (d)] in the 
spasing state. 
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FIG. 4: Ultrafast dynamics of spaser. (a) For monostable 
spaser (without a saturable absorber) , dependence of SP pop- 
ulation in the spasing mode N n on time t. The spaser is sta- 
tionary pumped at a rate of g = 5 x 10 12 s _1 . The color-coded 
curves correspond to the initial conditions with the different 
initial SP populations, as shown in the graphs, (b) The same 
as (a) but for the temporal behavior of the population inver- 
sion U2i- (c) Dynamics of a monostable spaser (no saturable 
absorber) with the pulse pumping described as the initial in- 
version ri2i = 0.65. Coherent SP population N n is displayed 
as a function of time t. Different initial populations are indi- 
cated by color-coded curves, (d) The same as (c) but for the 
corresponding population inversion 7221. (e) The same as (a) 
but for bistable spaser with the saturable absorber in concen- 
tration p a = 0.66p. (f) The same as (b) but for the bistable 
spaser. (g) The same as (e) but for the pulse pumping with 
the initial inversion 7221 = 0.65. (h) The same as (g) but for 
the corresponding population inversion 7221. 



IV. SPASER AS ULTRAFAST QUANTUM 
NANOAMPLIFIER 

A. Problem of Setting Spaser as an Amplifier 

Now we consider the central point of this article: set- 
ting the spaser as an ultrafast quantum nanoamplifier. 
As we have already mentioned in the Introduction (Sec. 
HJ, the principal and formidable problem is that, in con- 



trast to the conventional lasers and amplifiers in quan- 
tum electronics, the spaser has an inherent feedback that 
cannot be removed. Thus the spaser will always develop 
generation and accumulation of the macroscopic number 
of coherent SPs in the spasing mode. This leads to the 
the population inversion pinning in the CW regime at a 
very low level - cf. Fig. fj] (b). This CW regime corre- 
sponds to the net amplification equal zero, which means 
that the gain exactly compensates the loss, which condi- 
tion is expressed by Eq. (fl~2"|) . This is a consequence of the 
nonlinear gain saturation. This holds for any stable CW 
generator (including any spaser or laser) and precludes 
using them as amplifiers. 

We propose two regimes for setting the spaser as a 
quantum nanoamplifier. The first is a transient regime 
based on the fact that the establishment of the CW 
regime and the consequent inversion pinning and total 
gain vanishing requires some time that is determined 
mainly by the rate of the quantum feedback, but de- 
pends also on the relaxation rates of the SPs and the 
gain medium. After the population inversion is created 
by the onset of pumping and before the spasing sponta- 
neously develops, as we show below in this Section, there 
is a time interval of approximately 250 fs, during which 
the spaser provides usable (and as predicted, quite high) 
amplification. 

The second way to set the spaser as a quantum 
nanoamplifier is a bistable regime that is achieved by 
introducing a saturable absorber, which prevents the 
spontaneous spasing. Then injection of a certain above- 
threshold amount of SP quanta will saturate the absorber 
and initiate the spasing. Such a bistable quantum ampli- 
fier will be considered in the subsequent Subsection. The 
bistability is the most promising regime which allows not 
only the functioning of the spaser as a ultrafast threshold 
(logical) amplifier, but also a quasi-CW bistability where 
the spaser works as a memory element [cf. Sec. IIII B| 
with an ultrafast switching time. 

The temporal behavior of the spaser has been found by 
direct numerical solution of Eqs. ([I])-©. This solution 
is facilitated by the fact that in the model under consid- 
eration all the chromophores experience the same local 
field inside the nanoshell, and there are only two types of 
such chromophores: belonging to the gain medium and 
the saturable absorber, if it is present. 



B. Monostable Spaser as a Nanoamplifier in 
Transient Regime 

In this Subsection we consider a monostable spaser in 
a transiebt regime. This implies that no saturable ab- 
sorber is present. We will consider two pumping regimes: 
stationary and pulse. 

Starting with the stationary regime, we assume that 
the pumping at a rate (per one chromophore) of g = 
5 x 10 12 s _1 starts at a moment of time t = and stays 
constant after that. Immediately at t = 0, a certain num- 
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ber of SPs are injected into the spaser. We are interested 
in its temporal dynamics from this moment on. 

The dynamical behavior of the spaser under this pump- 
ing regime is illustrated in Figs. [4] (a), (b). As we see, 
the spaser, which starts from an arbitrary initial popula- 
tion N n , rather rapidly, within a few hundred femtosec- 
onds approaches the same stationary ( "logical" ) level. At 
this level, an SP population of N n = 67 is established, 
while the inversion is pinned at a low level of rii\ = 0.02. 
On the way to this stationary state, the spaser experi- 
ences relaxation oscillations in both the SP numbers and 
inversion, which have a trend to oscillate out of phase 
[compare panels (a) and (b)]. This temporal dynamics 
of the spaser is quite complicated and highly nonlinear 
(unharmonic) . It is controlled not by a single relaxation 
time but by a set of the relaxation rates. Clearly, among 
these are the energy transfer rate from the gain medium 
to the SPs and the relaxation rates of the SPs and the 
chromophores. 

In this mode, the main effect of the initial injection of 
the SPs (described theoretically as different initial val- 
ues of N n ) is in the interval of time it is required for 
the spaser to reach the final (CW) state. For very small 
N n , which in practice can be supplied by the noise of the 
spontaneous SP emission into the mode, this time is ap- 
proximately 250 fs (cf.: the corresponding SP relaxation 
time is less then 50 fs). In contrast, for the initial values 
of N n = 1 — 5, this time shortens to 150 fs. Whether 
this is a practically usable difference remains to be seen 
at more advanced stages of development. 

Now consider the second regime: pulse pumping. The 
gain-medium population of the spaser is inverted at t = 
to saturation with a short (much shorter than 100 fs) 
pump pulse. Simultaneously, at t = 0, some number of 
plasmons are injected (say, by an external nanoplasmonic 
circuitry). In response, the spaser should produce an 
amplified pulse of the SP excitation. Such a function of 
the spaser is illustrated in Figs. [4] (c) and (d). 

As we see from panel (c), independently from the initial 
number of SPs, the spaser always generates a series of SP 
pulses, of which only the first pulse is large (at or above 
the logical level of N n ~ 100). (An exception is a case 
of little practical importance when the initial N n = 120 
exceeds this logical level, when two large pulses are pro- 
duced.) The underlying mechanism of such a response 
is the rapid depletion of the inversion seen in panel (d), 
where energy is dissipated in the metal of the spaser. The 
characteristic duration of the SP pulse ~ 100 fs is defined 
by this depletion, controlled by the energy transfer and 
SP relaxation rates. This time is much shorter than the 
spontaneous decay time of the gain medium. This accel- 
eration is due to the stimulated emission of the SPs into 
the spasing mode (which can be called a "stimulated Pur- 
ccll effect"). There is also a pronounced trend: the lower 
is initial SP population N n , the later the spaser produces 
the amplified pulse. In a sense, this spaser functions as 
a pulse-amplitude to time-delay converter. 



C. Bistable Spaser (with Saturable Absorber) as 
an Ultrafast Nanoamplifier 

Now let us consider the most important point of this 
article: a bistable spaser as a quantum threshold ("log- 
ical") nanoamplifier. Such a spaser contains a saturable 
absorber mixed with the gain medium with parameters 
indicated at the end of Sec. Ill Al In particular, the con- 
centration of the saturable absorber p a — 0.66p. This 
case of a bistable spaser amplifier is of a particular in- 
terest because in this regime the spaser comes as close 
as possible in its functioning to the semiconductor-based 
(mostly, MOSFET-based) digital nanoamplifiers. As in 
the previous Subsection, we will consider two cases: the 
stationary and short-pulse pumping. 

We again start with the case of the stationary pumping 
at a rate of g = 5 X 10 12 s _1 . We show in Figs. 0] (e), 
(f) the dynamics of such a spaser. For a small initial 
population N n = 5 x 10~ 3 simulating the spontaneous 
noise, the spaser is rapidly (faster than in 50 fs) relaxing 
to the zero population [panel (e)] , while its gain- medium 
population is equally rapidly approaching a high level 
[panel (f)] ri2i = 0.65 that is defined by the competition 
of the pumping and the enhanced decay into the SP mode 
(the purple curves). This level is so high because the 
spasing SP mode population vanishes and the stimulated 
emission is absent. After reaching this stable state (which 
one can call, say, "logical zero"), the spaser stays in it 
indefinitely long despite the continuing pumping. 

In contrast, for initial values N n of the SP population 
large enough [for instance, for N n — 5, as shown by the 
blue curves in Figs. [4] (e) and (f)], the spaser tends to the 
"logical one" state where the stationary SP population 
reaches the value of N n « 60. Due to the relaxation 
oscillations, it actually exceeds this level within a short 
time of < 100 fs after the seeding with the initial SPs. As 
the SP population N n reaches its stationary (CW) level, 
the gain medium inversion 77,21 is pinned down at a low 
level of a few percent, as typical for the CW regime of the 
spaser. This "logical one" state salso persists indefinitely, 
as long as the inversion is supported by the pumping. 

There is a critical curve (separatrix) that divide the 
two stable dynamics types (leading to the logical levels 
of zero and one) . For the present set of parameters this 
separatrix starts with the initial population of N n ss 1. 
For a value of the initial N n slightly below 1 , the SP pop- 
ulation N n experiences a slow (hundreds fs in time) relax- 
ation oscillation but eventually relaxes to zero [Fig.|4](e), 
black curve], while the corresponding chromophore pop- 
ulation inversion 77.21 relaxes to the high value 77.21 = 0.65 
[panel (f), black curve]. In contrast, for a value of N n 
slightly higher than 1 [light blue curves in panels (e) and 
(f )] , the dynamics is initially close to the separaratrix but 
eventually the initial slow dynamics tends to the high SP 
population and low chromophore inversion through a se- 
ries of the relaxation oscillations. The dynamics close to 
the separatrix is characterized by a wide range of oscil- 
lation times due to its highly nonlinear character. The 
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initial dynamics is slowest (the "decision stage" of the 
spaser bistable that lasts > 1 ps). The "decision time" 
is diverging infinitesimally close to the separatrix, as is 
characteristic of any threshold (logical) amplifier. 

The gain (amplification coefficient) of the spaser as a 
threshold (logical or bistable) amplifier is the ratio of the 
high CW level to the threshold level of the SP population 
N n . For this specific spaser with the chosen set parame- 
ters, this gain is ~ 60, which is more than sufficient for 
the digital information processing. Thus this spaser can 
make a high-gain, ~ 10 THz-bandwidth logical ampli- 
fier or dynamical memory cell with excellent prospects of 
applications. 

The last but probably the most important regime to 
consider is that of the pulse pumping in the bistable 
spaser. In this case, the population inversion {n^i = 
0.65) is created by a short pulse at t = and simulta- 
neously initial SP population N n is created. Both are 
emulated as the initial conditions in Eqs. (UJ)-©. The 
corresponding results are displayed in Figs. [4] (g) and 
(h). 

When the initial SP population exceeds the critical 
one of N n = 1 (the blue, green, and red curves), the 
spaser responds with generating a short (duration less 
than 100 fs) pulse of the SP population (and the corre- 
sponding local fields) within a time < 100 fs [panel (g)]. 
Simultaneously, the inversion is rapidly (within ~ 100 fs) 
exhausted [panel (h)]. 

In contrast, when the initial SP population N n is less 
than the critical one (i.e., N n < 1 in this specific case), 
the spaser rapidly (within a time < 100 fs) relaxes as 
N n — > through a series of realaxation oscillations - see 
the black and magenta curves in Fig. 0] (g). The corre- 
sponding inversion decays in this case almost exponen- 
tially with a characteristic time ~ 1 ps determined by the 
enhanced energy transfer to the SP mode in the metal - 
see the corresponding curves in panel (h). Note that the 
SP population decays faster due to the threshold nature 
of spasing. 

Gold as the spaser nanoplasmonic-core metal, which 
is used in experiments^ 6 -, has the relaxation rate an or- 
der magnitude higher than silver. Therefore it is plau- 
sible that a gold-core spaser as an amplifier will have a 
bandwidth ^100 THz. We will consider such a system 
elsewhere. Taking into account that the local fields pro- 
duced by the spaser are concentrated within a few tens of 
nanometers in its vicinity, these are parameters that bear 
promise of many applications in the fundamental science 
and coming digital and analog femtosecond technologies. 



V. DISCUSSION AND CONCLUSION 

We have demonstrated principal possibilities of the 
spaser functioning as a ultrafast nanoamplifier of the lo- 
cal optical fields. In doing so we have had to overcome 
a principal and formidable problem. Any spaser has the 
inherent feedback. Consequently, the spaser as proposed 



initially^ and observed experimentally 16 rapidly, on the 
femtosecond time scale, will generate the coherent popu- 
lation of the SPs and approach the CW regime. In this 
regime, the net gain (the saturated amplification minus 
loss) is exactly zero (the same as for the lasers). Thus CW 
spasers (or CW lasers, for that matter) principally cannot 
be used as amplifiers. In the conventional optical ampli- 
fiers, in contrast to lasers, the feedback is deliberately 
removed and any parasitic feedback (scattering, etc.) is 
carefully minimized to prevent the spontaneous genera- 
tion. Such a regime is not possible for the spaser as it is 
presently known. In this article, we have demonstrated 
two ways to set a spaser as a nanoamplifier: (i) the con- 
ventional monostable spaser can amplify in the transient 
ultrafast (femtosecond) regime before the CW generation 
is established; (ii) the spaser with the saturable absorber 
can become bistable and function as threshold (logical) 
amplifier with ultrafast (femtosecond) switching times. 

We have developed quantum theory of the spaser based 
on density matrix (optical Bloch) equations. The spaser 
described is a nanoparticle that consists of the metal core 
and the gain (active) medium. The metal core plays the 
role of the resonator ( "cavity" ) whose states are the local- 
ized SPs, which we have quantized. The gain medium is 
described completely quantum mechanically taking into 
account the excitation, decay, and dephasing processes. 
The main approximation that we have employed is that 
the SP population number N n is treated as a classical 
variable. This is the same semiclassical approximation 
that is most often used in the physics of lasers. Given 
the fact that the typical number of the SP quanta per 
spasing mode is still large enough (N n ~ 100), this ap- 
proximation appears reasonable and reliable. 

Pursuing the goal of the fastest possible operation and 
the nanometric size (on the same order as that of micro- 
electronic MOSFETs), we have deliberately considered 
spasers whose size is much less than the radiation wave- 
length and whose metal thickness is less than the skin 
depth. This has allowed us to use the quasistatic SP 
eigenmodes as full counterparts of the electrostatic fields 
in the transistors. Note that such nanoscopic a spaser 
has recently been demonstrated experimentally^ 6 -. 

Though the main emphasis of this article is on the ul- 
trafast dynamics of the spaser as a nanoamplifier of the 
local plasmonic fields, the CW regime is described first 
as the basis of the understanding the ultrafast dynamics. 
The corresponding results are illustrated in Fig. [5] The 
"spasing curve" (the dependence of the SP population 
N n on the pumping rate) shows a pronounced threshold 
and a linear dependence of the coherent SP population 
N n on the pumping rate g [Fig. [2] (a)]. This linear depen- 
dence N n (g) is due to the extremely strong feedback in 
the spaser caused by the very strong SP mode localiza- 
tion. The same effect is responsible for clipping the pop- 
ulation inversion in the spasing state to very low levels, 
as illustrated in Fig. [2] (b) . With the increased pumping, 
a higher SP population leads to lower quantum fluctu- 
ations and decreased spectral width of the spasing line, 
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as shown in Fig. [5] (c). Gradually, a very narrow and in- 
tense spasing line appears in the emission spetrum of the 
spaser - see Figs. [2] (d)-(f). This narrow line is indica- 
tive of the intense optical fields that are excited in the 
vicinity of a spaser - cf. Fig. QJa) where the field must 
be multiplied by a factor of \/N n . Qualitatively, the lin- 
ear spasing curve and the narrowing of the spasing line 
with the pumping are in a full agreement with the recent 
experiment^. 

Adding a saturable absorber to the spaser shell, we 
have shown that it can work as a bistable based on quan- 
tum amplification. The corresponding results are illus- 
trated in Fig. [3] for a CW operation mode. 

The main results of this article concern with ultrafast 
processes of amplification in the spaser. These are illus- 
trated in Fig. 01 We have shown that a monostable spaser 
(without a saturable absorber) performs amplification of 
the injected SPs in a transient regime. However, this 
process only manifests itself in the shorter times to es- 
tablish the high SP population - see Figs. 0] (a)-(d) and 
the corresponding discussion. In this sense, such a spaser 
performs the function of the amplitude to time converter 
but not a true nanoamplifier. 

The most significant results are those related to the 
ultrafast amplification in the bistable spaser containing 
a saturable absorber within its gain shell. This ampli- 
fication has many similarities with that in semiconduc- 
tor threshold (logical) nanoamplifiers based on the MOS- 
FETs. Such a spaser with a stationary pumping behaves 
as an ultrafast nanoscale bistable - see results in Figs. 
0] (e) and (f). For the initial conditions not too close to 
the threshold, it amplifies and switches with a character- 
istic time of ~ 100 fs. The bistable spaser with the pulse 
pumping is a high-gain ultrafast pulse nanoamplifier - 
see Figs. 0] (g) and (h) and the corresponding discussion 



in the previous Section. Its typical response time and 
pulse length are < 100 fs, which correspond to the 10 
THz bandwidth. 

For the gold-based spaser, in contrast to the present 
silver-based one, the typical switching time should be 
much shorter, presumably ~ 10 fs corresponding to ~ 100 
THz bandwidth. However, the gold-based spasers will 
work at a lower SP population due to the higher loss in 
the metal. Correspondingly, they will possess lower sig- 
nal/noise ratio and wider spectral lines. We will consider 
such spasers elsewhere. 

Concluding, the spaser is the size of the MOSFET 13 
(~ 10 nm) and can perform the same function of signal 
amplification on the nanoscale. Having a high enough 
gain, the spaser can be an active-element foundation of 
highly integrated nanoplasmonic devices, including ultra- 
fast processors. Based on metals, it is inherently much 
faster, by a factor of ~ 1000 than silicon devices, as we 
have shown above. For the same reason, it also highly ro- 
bust environmentally: it can work at high temperatures, 
in the presence of microwave and ionizing radiations, etc. 
This article has presented quantum theory of the spaser 
as a nanoamplifier, which is critically needed for under- 
standing and utilization of the spaser. 
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